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Abstract 

The Casimir force between parallel plates of arbitrary kind is shown 
to be simply related to the plates transmission and reflection coeffi- 
cient. A trivial application of this general relation leads to known 
Lifshitz force between dielectrics as well as its generalizations. 



The simple two ideally conducting parallel plates Casimir force F 
"5ioc?0 recently been measured experimentally @, [10| . It has various gen- 



eralizations to different kinds of fields, boundary conditions and geometries 
whose calculation is in general quite complicated^, [7], §@. A new vari- 
ant in which each of the plates conducts in specific in general different di- 
rection has been exactly evaluated recently^] and can hopefully be seen 
experimentally [0]. Another generalization is to replace the ideal conducting 
plates by dielectrics. Thus we may for example consider the vacuum energy 
corresponding to electromagnetic action in matter / d 4 x(eE 2 — j^B 2 ) as 
some very complicated functional of e(x),fi(x). In general this energy will 
contain divergent terms. When we consider the energy difference between 
two configurations (given by e(x),fi(x)) the divergent terms should cancel 
provided that the configurations may be related by a physical process which 
does not involve a change in internal matter properties (otherwise such a 
change would bring about a dependence on the physical UV cutoff) . In par- 
ticular relative displacement of the different dielectrics in space should result 
in a (finite) force. The special case of the force between two infinite parallel 

f ei z < 

dielectrics (corresponding to e = < 1 < z < a ) was considered long ago 

[ e 2 a < z 

by Lifshitz||. 

We shall consider in the following general configurations depending only 
on the z coordinate. We shall show using path integral techniques that in this 
case the vacuum Casimir energy may be simply related to the transmission 
coefficient through the system This reduces the computation of the force to a 
solution of an elementary one dimensional scattering problem. In particular 
to find the electromagnetic Casimir force one merely need to know the clas- 
sical reflection and transmission coefficient through the relevant system. In 
fact the relation between the Casimir energy and the transmission coefficient 
follows from of the relation between scattering amplitudes singularities and 
the eigenenergies. 

We use the functional integral formalism to write the Casimir energy as 
e~ ET = J T><fte~ s = det(...)~3. Consider first a scalar field described by an 
action S — j <i 4 x(0D0 + e(z)<f) 2 ). To compute the functional determinant 
of this we first Wick rotate and Fourier transform in the three directions 
in which there is a translation symmetry-t, x, y. After doing that we are 
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left with the action / d 3 k J dz{—(j)d 2 z (j) + (e(z) + k 2 )4> 2 ). The total Casimir 
energy for this problem is given by summing contributions corresponding to 
all possible values of k = (k t , k x , k y ). 

E = \A J log det Hk where detTCk is the functional determinant of 
the one dimensional schrodinger like operator Ti^ = —d 2 + e + k 2 . 

Discretizing this using J dz<p(k 2 -d 2 z +V(z))<p ~ £ ((^ 2 + V(-Zi))^ - / I+1+ ^ 2 ' 2 1 " 20 ' ) Az 
leads us to the determinant 

(k 2 + V( Zl ))Az 2 + 2 -1 . . . 

-1 (fc 2 + ^(^ 2 ))Az 2 + 2 -1 

-1 



. . . (k 2 + V{z n ))Az 2 + 2 

The determinant easily /„ is seen to satisfy a recursion relation I n = ((k 2 + 
V(z n ))Az 2 + 2J n _! - J n _ 2 

which may also be written in the form Jn+Jn ^ 2 2 ^ 2Jn ~ 1 = (k 2 + V(z n ))I n _i 

hence in the continuum limit I" — (k 2 + V)I = 

i.e. the determinant satisfies the corresponding schrodinger equation. 

This is in complete accordance with another well known trick for cal- 
culating determinants of one dimensional schrodinger operators. 

We are assuming that e(z) vanishes at z — ±oo. Under this assumption 
we know that the solutions of Titp = behave asymptotically as ip ~ e ±kz . 
Let us look for a solution which is regular at minus infinity i.e. ip ~^ ekz f° r 
z — > — oo. Such a solution always exist and is unique .However in general it 
will not be regular at +co. Rather it behaves in this limit as if) — > Ae kz + 
Be~ kz . The function if} is a regular solution of TCif> = if and only if the 
constant A(k) = limit z ^ OQ (ip(z)e~ kz ) is zero. On the other hand we know 
that a regular solution (which is necessarily given by the above defined ifi) 
exists if and only if detTi = 0. This suggests the identification detH = 
\im(ifj(z)e kz ). 

A more rigorous treatment starts from a similar problem where the parti- 
cle is confined to a finite segment [zi, z^\ of length L with boundary conditions 
ip(zi) = ip(z2) = 0. In this case we may define a function ip(z) to be the 
solution of TCif> = satisfying the boundary conditions if){zi) = Q,if)'(zi) = 1 
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and we easily see that ip(z2) = <^4> det Tt — 0. The advantage of dealing first 
with the problem on a finite segment is that now the argument is correct for 
complex valued k while the previous argument was justified only if Re(k) > 
(while we know that most eigenvalues of 7i have in fact Re(k) = 0). Hence 
we may now conclude that ipk{ z 2) and det Tik have exactly the same zeroes all 
over the complex plane. We also know from general properties of schrodinger 
operator in one dimension that all the zeroes of det Ti are simple since all 
eigenvalues are nondegenerated.As k — >• oo the problem becomes classical 
and one should be able to deduce the behaviour of both ^(^2) and of det Tik 
at that limit from semi classical considerations. This allows concluding that 
ipk(z2) and det TCk are equal up to an (e(a;)-independent) normalization con- 
stant. 

We are interested in knowing det H only for real positive k. For such k 
the boundary conditions ip(zi) = may be safely replaced by demanding 
ip — > e kz for z — > —00. Making this change just amounts to neglecting a term 
of the order e~ kL which disappears once when L — >• 00. By the same reason 
we may also replace ip{z2) by the coefficient of the increasing exponential e kz 
in the asymptotic expression for ip(z) (at this step we also renormalized the 
determinant by the constant factor e kL which is legitimate). It is also easy 
to see that had we used different boundary conditions such as ip'{zi) = 
instead of ip(zi) = the result would only change by something proportional 
to e~ kL . 

Given any function e(z) we may according to the above compute det 7"^ 
by solving (possibly numerically) a one dimensional schrodinger equation and 
then find the Casimir energy from the relation E = J log det TLk- We 
saw that det Tik is given by the coefficient of e kz in the expression of ip at 
the right of the potential(i.e. for z 3> 1) when we assume that left of the 
potential (i.e. z <^ —1) ip was given by purely e kz . Analytically continuing 
this to imaginary values of k writing k = —iq we see that the analytically 
continued detHk becomes just the (inverse) transmission coefficient for a 
wave of momentum -q hitting the potential e(z) from the right to propagate 
into the other the side of the potential. This is of course identical to the 
more conventionally defined (inverse) transmission coefficient for a wave of 
momentum +q coming from the left to propagate into the right of the poten- 
tial. Denoting by t(q) the transmission coefficient and changing back from q 
to k we have the result: detHk = \t(ik)Y l . 
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Next we would like to study the case where the potential e(z) is generated 
by two objects one located around z = and the other around z = a thus 
obtaining their mutual force. 

Scattering in one dimension is usually described by a 2 x 2 transfer matrix 

! i \ (satisfying tt* + rr* = A wave of the form 

Ae iqz + Be~ iqz coming from the left will be transformed upon crossing the 

scattering potential to the right into Ce iqz +De~ iqz where ^£)^ = ^^^^- 

As the location of the scattering potential is moved by a distance a to the 
right without changing its shape the transfer matrix T change into A~ 1 TA 

( e iqa \ 

where A is the matrix representing translation: A = I _ iqa I . 

For a potential e(z) generated by two different slabs one located around 
z = and the other around z = a the total transfer matrix of the system 
is given by the product of the two transfer matrices corresponding to the 
two objects: T\ and A~ 1 T 2 A. so that we have T = A~ l T 2 ATi. And in 
particular r 1 = (T) 2 , 2 = (T 1 ) 1)2 (T 2 ) 2il e 2i ^ + (T 1 ) 2 , 2 (T 2 ) 2 , 2 = |ge 2 ^ + ^. 

Assuming real potential we have the relations r(g)* = r(—q),t(q)* = t(—q) 
Which allow us to make the analytic continuation into q = ik obtaining 
det Hk = ^(lifc)^^) e~ 2fca + tl ^ t2 ^ ■ For the Casimir force of interest acting 
between the bodies we need only the dependence on the distance a. Thus 
we may renormalize defHk into 1 + ri ~^~r~^fc e~ 2ka . Hence the Casimir 
interaction between the two bodies will be given by 

1 f d 3 k i ri(—ik)r 2 {ik)ti(ik) _ 2k 



E =2iWf l ° i(1+ ti(-ik) « 



For a massive field one obtains similar relation with e~ 2ak — > e ~ 2a ^ k ' 2 + m ' 2 . 
We remark that upon replacing a potential by its mirror image e(—z) 
the reflection coefficient r changes into f = —rt*t~ l hence the last result may 
be also written as 

1 r rl 3 k 

E= 2j^f l0g(1 " ri W) e '*) ( 2 ) 



^^The ordinary S-matrix is of the form S = ( „ 
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which is explicitly symmetric under T\ ^> T 2 _1 as expected. f\ 

As a simple example suppose that the scattering potential are given by 
delta-functions: e = \i5(z)+\25(z— a). The transfer and reflection coefficient 
from a delta potential are t(q) = 2 ^ A , r(q) = 2i _ y Substituting this in the 

general result(||) we have E = | / log(l — Ai ^ 2fc x^+2k e ~ 2ka ^- ^ n P ar ticular 
we see that for A — > oo this reduces as expected to the standard Casimir 
energy of two parallel plates | J log(l — e ~ 2ka ) = — 14 ^ a3 ■ We can also 
see that when A becomes negative the energy develops an imaginary part as 
expected. It is also quite easy to generalize the last result to a situation where 
there are more then two plates Thus obtaining an explicit expression for a 
three body force which should exist as long as there is no perfect reflection 
from the middle plate. Also the case of several finite width slabs including 
circular periodic boundary conditions can be addressed by these methods O 



The fact that the Casimir force may be expressed in terms of reflection 
and transmission coefficients seems very natural. It is known that one way for 
obtaining the Casimir force is by summing the pressure exerted on the various 



bodies as a result of the vacuum zero modes being scattered of them [^2 
Clearly such an approach should also result in expressing the force in terms 
of reflection and transmission coefficients. 

If instead of a single scalar we have a multicomponent field most of 
the above arguments can be repeated. The only difference being that now 
t(k),r(k) are matrices in internal space and that we have to take the (fi- 
nite dimensional) determinant also with respect to them. Thus E/A = 
— | / log det t(ik) in particular if the different components of the field 
do not mix the determinant becomes just the product of the various factors 
and the energy the sum of the different energies as expected. 

In order to apply this to the Casimir force between two parallel di- 
electrics of dielectric constants £1,62 and widths Di,D 2 we just need to 
solve the simple classical problem of computing the transmission of a plane 
electromagnetic wave through the system made of these two plates. The 
classical transfer matrix for an electromagnetic wave when passing from 
one dielectric to another through the z = plane can be expressed quite 

2 If the two slabs are mirror image of eachother i.e 7*1 = ri then eq(||) imply that they 
would always attract 

3 Also for this case we can consider arrays of parallel slabs. For the case mentioned 
above of plates conducting in different directions which is relevant only for vectorial fields 
there is partial reflection and incomplete screening of nonadjacent plates 
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simply in terms of the ratio ei i2 = ei/e 2 and the ratio «i j2 = Qz /Qz 



For a wave whose polarization is parallel to the z = plane 




we have = J-( a ^ 01 \ \ while for the other polarization = 
2a \a-la + lJ y 

\ ( e ix \ 

I . Let us also denote as before A(x) = I ^ ^_ ix I . 

If we have a system of two parallel dielectrics of widths Di, D 2 and dielec- 
tric constants ei,e 2 separated by a distance a inside a medium of dielectric 
constant €q then the total transfer matrix through the system is given by 

T = A{-q z {D x + D 2 + a))T(e , e 2 )A(q^ D 2 )T(e 2 , e )A(q z a)T(e , e 1 )A(g( 1 ) J D 1 )T(e 1 , 

(3) 

The contribution of each polarization to the Casimir force between the two 
plates is | J log(T 2i2 (g z = zA;)). Since the two polarization are completly 
independent the total force is just the sum of these two contributions. It 
is more convinient to make the changes u — > ik t and q z —> ik already in 
the expressions ©for the matrix T This amounts to the changes A(x) — > 

( e e x ) ,a ^ and p° ssibl y also e M e (^*)- 

Substituting all this into eq(|3|) one easily finds that the correspond- 
ing functional determinant for the two polarizations (up to unimportant 
a-independent factors) are: 



ei + ai + 2e 1 aicoth(a 1 D 1 ) jy el + a% + 2e 2 a 2 coth(a 2 D 2 ) J 

(4) 



det ^ e ) = 1 - ( 1 ~ a ^ )( l ~ a i ) e 

k 1 + a{ + 2a 1 coth(a 1 D 1 ) JK l + a 2 + 2a 2 coth(a 2 D 2 ) 



-2ka 



(here we assumed e = 1 and hence a, = 

The plate interaction energy is given (up to a-independent terms) by E = 
i/(0(logdet^ 0) + logdet^ e) ). 

In particular if the widths Di,D 2 are large (compared to a) then (£|) re- 
duces to det^i 0) = l-(^) 0jl (^) 0i2e - 2fca and detH^ = l-(^)o,i(S)o, 2 e- 2fca 
for the two polarizations. Substituting this into the general formula leads to 
the well known Lifshitz force H. 
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Equation (^) may also be used in more general scenarios e.g. to find 
exactly the torqe between two nonisotropic dielectric plates (evaluated in 
some approximation in||ll||). The resulting expression for the energy may 



turn to be quite complicated in such a general case. Some simplification 
always occur for large widths Di,D 2 since than A(q z D) reduces in@ to 

the projection ^ ^ therefore ignoring a-independent factors we have 

The above results reduce the problem of computing the Casimir energy 
for general plate geometry to the problem of computing the reflection coef- 
ficient from the corresponding one dimensional potential. Generalizing this 
method to arbitrary non planar geometry would be very difficult if not im- 
possible. However one may still try to use similar method for problems of 
spherical symmetry. Using partial wave expansion reduces this to an essen- 
tially one dimensional problem. Thus we may identify (each specific partial 
wave contribution to) the functional determinant with the ratio of the out- 
going to the ingoing amplitude i.e. with the phase shift or more accurately 
with e 2lSl . This consideration leads us immediately to the following exact 
formal expression for the Casimir energy 

£ = E(2/+i)/ o (5) 

This result also follows imeadiatly from known relation between phase shifts 
and density of states n(E) = --f§- The same result can also be obtained 
by more standard considerations. Suppose we want to compute the Casimir 
energy by explicitly finding the eigenfrequncies and summing | Tw- Since 
the asymptotic form of the partial wave at infinity is ~ - sin(/cr — df) we 
easily see that enforcing a boundary condition stating that the fields should 
vanish on some sphere of very large radius R leads to the simple equation 
kR — Si(k) = nir for the eigenvalues. The Casimir energy of a specific partial 
wave is therefore Ei — ~ J2n k = \ nn+ ^ l( - kS) = const + In this 

equation k = k(n) is the solution of kR — 5i(k) = nir however since the di- 
ameter a of the scattering potential is much smaller then R we always have 
kR » 5i(k) and therefore we may write simply k — As R » a the 
sum over n may be replaced by an integral J2n ~^ R Jo°° ^ lading us back 
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to eq(^|). However equation^) unlike(^) is in general divergent and needs a 
nontrivial regularization. The difficulty originates from the fact that our pre- 
vious normalization condition that of demanding the energy to vanish as the 
different bodies are moved infinitly apart from eachother without changing 
their shape becomes meaningless here. 



An issue related to the above which we elaborate elsewhere [jig] is the 
question of experimental testing the repulsive nature of the Casimir energy 
of the sphere || |7|]. One cannot contrary to what has been (implied/)stated by 
various authoursfl3| || [To]] check this by finding that two near by conducting 



hemispheres repel. Indeed the latter case of saperated hemispheres represents 
a well defined and in principle testable experimental set-up. However such 
experiment are likely to manifest attraction of the two hemispheres. Thus 
consider the two thin rims at the boundaries of the hemispheres which are 
a small distance 5 apart. On their own the rims have an attractive Casimir 
attraction of F rims ~ 2'kRM. Since the wavelength of vacuum fluctuations 
which generate most of this attraction are A ~ 5, the existance of the remain- 
ing bulk of the hemispheres of radius R 5 is not expected to drastically 
modify this force. Since this force is vastly larger (by (R/5) 3 \) than the 
repulsive force of a single sphere F sp h er e ~ |§ the attractive nature is not 
modified. 

Acknowledgment I thank Prof S.Nussinov for discussions and his in- 
terest. 
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